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We investigate the possibility of searching for ultra high energy neutrinos in cosmic rays using 
acoustic techniques in ocean water. The type of information provided by the acoustic detection is 
complementary to that of other techniques, and the filtering effect of the atmosphere, imposed by the 
fact that detection only happens if a shower fully develops in water, would provide a clear neutrino 
identification. We find that it may be possible to implement this technique with very limited 
resources using existing high frequency underwater hydrophone arrays. We review the expected 
acoustic signals produced by neutrino-induced showers in water and develop an optimal filtering 
algorithm able to suppress statistical noise. The algorithm found is computationally appropriate to 
be used as a trigger for the signal processors available on existing arrays. We estimate the noise 
rates for a trigger system on a very large size hydrophone array of the US Navy and find that, while 
a higher density of hydrophones would be desirable, the existing system may already provide useful 
data. 

PACS numbers; 95.85.Ry, 13.85.Tp, 96.40.Tv, 96.40.-z 



I. INTRODUCTION 

The understanding of the highest energy cosmic rays 
represents one of the most challenging fields of modern 
physics. While abundant high quality data has greatly 
enhanced our knowledge of cosmic rays of energies up to 
several tens of GeV, the study of higher energies is lim- 
ited by the low fluxes available. Yet the study of parti- 
cles with energies in excess of 10"'^^ eV (ultra-high energy 
cosmic rays, or UHECR) promises boundless opportu- 
nities of discovery. To date some 12 cosmic ray events 
have been observed with energy in excess of lO^o eV [|. 
While the acceleration mechanisms at these energies are 
not completely understood even more fundamental 
problems arise from the apparent inconsistency of data 
with the Greisen-Zatsepin-Kuzmin (GZK) cutoff. Such 
cutoff ^ is expected to limit the maximum energy of 
protons of cosmological origin somewhere below 10^° eV, 
because of the finite (« 50 Mpc) inelastic collision length 
of such particles in the cosmic microwave background ra- 
diation (CMBR). Indeed the inelastic scattering of UHE 
protons off CMBR should at the same time suppress the 
proton flux and breed neutrinos from decay-products of 
pions. While the present data seem to indicate that 
photons and neutrinos are not the main component of 
UHECR reliable identification of the primary parti- 
cle type and its energy are essential parameters for the 
study of this problem. 

It was recently pointed out ||^ that a substantial 
UHE neutrino component may accompany UHE protons 
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and nuclei, due to neutrino production in cosmic beam 
dumps. Neutrino production at ultra-high energies ap- 
pears in fireball models of gamma-ray bursts [|[ §1 1 £^c- 
tive galactic nuclei ||, [To[ ^ |l3|, and, to a lesser 
degree, in Galactic mechanisms |15| and, as mentioned, 
the GZK process Another hypothetical source is 

the decay of heavy objects predicted by some theories, 
known as the "top-down" models of UHECR produc- 
tion |, 0, [l . 

Weakly interacting neutrinos could, unlike UHE 
gamma rays and protons, reach us from distant and 
powerful sources, opening a deeper horizon for astro- 
physics, cosmology and, possibly, high-energy particle 
physics 1^, ^ . While atmospheric neutrinos represent an 
irreducible background for Earth-based detectors, such 
background is expected to be modest because of the ex- 
tremely long decay length of pions at the energies of in- 
terest. 

Three techniques have been used until now to detect 
UHE cosmic-rays, all involving the showering of the pri- 
mary particle in the Earth's atmosphere. The shower 
is then detected either by observing the fluorescence or 
Cerenkov light induced by the ionizing tracks in the 
air, or by directly detecting the charged particles in 
the shower tail with scintillation counters scattered on 
the ground |o[ Ip. While some of the these tech- 
niques provide the largest acceptances obtained in par- 
ticle detectors, in general the study of UHECR is still 
hampered by the very low flux one has to be sensitive 
to. Typical fluxes are 
~ 1 km~^y~-'^ above 10^^ eV, and 
above 10^° eV Dedicated neutrino telescopes using 
Cerenkov light under water and the Antarctic ice-cap are 
for the time being optimized for the TeV to EeV energy- 
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region 

While much has been learned from the above detec- 
tors, it is important to explore, in parallel, new meth- 
ods that could either increase the flux sensitivity, and 
hence raise the energy threshold for detection, or help 
constraining the primary particle identification and its 
energy. Alternative methods being discussed include 
the detection of radio Cerenkov emission from the lu- 
nar soil (sensitive only to neutrinos that can cross 
the moon and interact upon exiting the satellite from 
its near side) and active radar detection of showers in 
the atmosphere (which has different systematics for 
particle identification since it is more sensitive to hori- 
zontal showers than other techniques). As discussed in 
this paper, another possibility consists in the detection of 
showers by means of the acoustic energy released in the 
medium where they develop. While, at least in principle, 
this technique could use as a radiator either the soil of 
the moon or a large body of water on the earth, here 
we concentrate on this second case that is more practi- 
cal and, as we will show, may be possible to test on a 
large scale with a very modest effort. It is important to 
realize that these two radiators would provide very dif- 
ferent information, the first being sensitive to UHECR 
of any type, while the second being sensitive essentially 
only to neutrinos because of the filtering effect of the 
earth's atmosphere. While the idea of taking advan- 
tage of the very high mechanical Q for small oscillation 
amplitudes and ultra-low seismic noise characteristic of 
the moon was first mentioned in , acoustic detection 
of ionizing particles in water was proposed in |Q and 
then developed |2q] in connection with the DUMAND 
project |2^, |2^, ^9[ . Although the primary goal of DU- 
MAND was optical Cerenkov detection of muon tracks in 
deep ocean water, the hydrophones, originally conceived 
to monitor photomultiplier positions, were proposed to 
be used for acoustic detection of neutrinos of > 10^^ eV 
energy. These early ideas were nicely complemented by 
experimental data collected at accelerators |30 . We will 
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review the past work on the subject and present a model 
useful to study UHE neutrino detection in sea water. We 
will then use this model to simulate signals and develop 
an algorithm to optimally filter hydrophone data to ex- 
tract the expected signals from statistical noise. While 
most the assumptions for the study are rather general, 
we will apply our results to the case of a large, high fre- 
quency test array that the US Navy operates off the coast 
of Florida. 

We note here that UHE neutrino fluxes substantially 
larger than the UHECR flux quoted by the experiments 
utilizing air showers could have gone unobserved. The 
neutrino cross section due to neutral and charged current 
interactions in the UHE regime is calculated in . Even 
for El, ~ 10^° eV the probability of interaction in travers- 
ing entire atmosphere's depth is only of the order of 
10-5. Many models |, |l^, ^ g predict 

neutrino fluxes higher than the flux of the UHECR ob- 
served in the atmosphere. For example, fireball gamma- 
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FIG. 1: Sound attenuation coefficient at 25°C, in the sea and 
distilled water as function of frequency |32|. 
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FIG. 2: Water temperature profiles for tlie tropical region of 
of the AUTEC array discussed below [Q. The two curves 
correspond to the extreme variations of temperature through 
the year. 



ray burst models predict a neutrino flux at energies 
El, > 10^° eV as high as ~1 km-^y"^, active galactic 
nucleus models 0.1 km-^y"^, and top-down models 
-10 km-^y-^, as quoted in fsl. 



II. ACOUSTIC SIGNATURE OF PARTICLES IN 
WATER 

At the energies of interest here, neutrinos interact in 
water by a deep inelastic scattering on quarks inside oxy- 
gen and hydrogen nuclei. The scattering produces a lep- 
ton and an hadronic shower, with similar energies shared 
among these two components. The structure of hadronic 
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showers is the same for all three flavors of (anti)neutrino. 
The behavior of the leptons, however, is different. In 
the case of t'c's the lepton energy goes in an electromag- 
netic shower and is essentially detected together with the 
hadronic energy. For i^^'s calculations by Mitsui ||3^ 
show that the mean-free-path in water between catas- 
trophic bremsstrahlung and direct pair production events 
with energy transfer greater than 10^^ eV is large (sev- 
eral km) with respect to the vertical size of the detectors 
considered here. The residual ionization along the track 
is low, so that muons are virtually undetectable by acous- 
tic methods Tau neutrinos create r leptons, whose 
mean free path turns out to be long enough to leave the 
detector volume for i?^^ > 10^^ eV. 

At El, ~ 10^" eV in water a hadronic shower deposits 
90% of its energy in a cylinder of some 20 cm radius 
and 20 m length. Sound is produced in water mainly 
by heating localized along the shower, resulting in vol- 
ume expansion, as first suggested in |2^ . The shower 
development (and hence the energy deposition) occurs 
at the velocity of light and can be regarded as instanta- 
neous for the purpose of acoustic phenomena. As already 
mentioned the sound generation has been confirmed ex- 
perimentally using artificial particles |Q . Although it is 
not completely clear to what extent these accelerator ex- 
periments, simulating '--^lO^" eV energies using bunches 
of over 10^^ protons of ~200 MeV, can describe the de- 
tails of the lower energy-density UHE neutrino interac- 
tions, their accuracy is probably sufficient for the present 
study. 

We analyze the acoustic signal production follow- 
ing ||2^. Let the energy deposited per unit volume per 
unit time be given by a function E{r,t). The total neu- 
trino energy is Eq = Jy E{r) d^r. The wave equation for 
the pressure pulse produced p is: 
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Where we use the parameters for sea water: c « 1500 m/s 
is the speed of sound, /3 « 1.2 x 10^"^ K^^ is the bulk co- 
efficient of thermal expansion, Cp w 3.8 x lO'^ J kg~^K"^ 
is the specific heat at constant pressure, and loq w 
2.5 X 10^" s~^ is the characteristic attenuation frequency. 
UJo is, strictly speaking, a function of frequency [Q, as 
can be seen from the plot of the attenuation coefficient 
a(dB/km) = (10^/ In 10woc)(27r/)2 in Figure |l|. For sim- 
plicity of calculations, we assume that coq is a constant 
in the frequency range / = 10 — 100 kHz characteristic 
of the signal. 

We note here that the coefficient of thermal expansion 
(3 depends upon the water temperature. It vanishes at 
~— 3°C for typical sea water of 3.5% salinity |3^. In the 
case of vanishing /3 other mechanisms of energy coupling 
to acoustic modes have been proposed |3^ . Extreme 
temperature profiles as a function of depth are given in 
Figure ^ for the tropical waters of the site discussed be- 
low. As we can see, although jS decreases with depth, it 
does not reach zero, so we will restrict ourselves to the 



case of thermal emission mechanism. 

The instantaneous nature of the heating mechanism 
can be expressed by E{r',t) = E{r')S{t). The pressure 
wave can then be calculated at location r as a function 
of time t as 

p(r,i) = / E{r')G{r-r',t)d^r' 
Jv 

where G(r, t) is the pressure pulse generated by a point 
source E = S{r)d{t), taking attenuation into account: 



where t ~ ^/r/JujQc). 

The region of the energy deposition is elongated in the 
direction of the initial velocity of the neutrino, which con- 
stitutes the axis of the cascade. The acoustic emission is 
coherent in the plane perpendicular to this axis. Thus, 
the radiation diagram has a pancake shape, perpendicu- 
lar to the shower axis, as shown in Figure |. This eflFect 
was also observed in the 30 cm long and 4.5 cm wide 
cylinder energy deposition by particles in one of the ac- 
celerator experiments described in | |30[ |. 

We calculate the time dependence of the pressure at a 
distance r — 1 km from the origin, for different directions 
of observation. As shown in Figure the angle 9 is calcu- 
lated in the forward direction from the perpendicular to 
the shower axis and the origin is at the starting point of 
the shower. Since the maximum energy deposition does 
not occur at the beginning of the shower, the maximum 
acoustic pulse results at some angle 9 > 0. 

To calculate the effect of a 10^° eV hadronic shower, 
we use a model based on the results of one-dimensional 
Monte Carlo simulations in ice . This model includes 
the interactions of 7r° and other short-lived resonances, 
as well as the Landau-Pomeranchuk-Migdal (LPM) ef- 
fect H, |o[ H, H, that are important at UHE. The 
transverse structure is modeled after . The results of 
acoustic emission by a hadronic shower are presented in 
Figure |4[ 

The pressure pulse has a typical bipolar shape and is 
shown, at the maximum of the radiation pattern, in Fig- 
ure |a. The peak pressure and energy fluence for other 
angles of observation are shown in Figure |^b. The signal 
at other observation angles has a smaller amplitude and 
is stretched in time, still preserving the bipolar shape. 

To calculate the case of an electromagnetic shower we 
use the Monte Carlo model LPMSHOWER [|| that in- 
cludes the LPM effect. In this case, with all the energy in 
the electromagnetic channel, the LPM effect dominates 
the shower shape, resulting in very large fluctuations, 
and a significant elongation for initial energies >1 EeV. 
For Ee = lO^'' eV the shower is ~300 m long. The 
non-uniform energy deposition is the result of individ- 
ual sub-showers starting at the locations of large energy 
loss. These sub-showers also create peaks in the acousti- 
cal radiation pattern. The properties of a typical shower 
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FIG. 3: Geometrical configuration of a shower and acoustic pulse produced. The angle 6 is calculated with respect to the 
forward direction from the perpendicular to the shower axis. The origin is at the starting point of the shower. 
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FIG. 4: Results of calculations of the acoustic signal from the hadronic part of the neutrino- induced shower ||3^, at the distance 
of 1000 m from the shower axis, for a primary hadronic energy of lO'^" eV: (a) the pulse shape at the observation angle of 
9 — 0.6°, where the amplitude is maximal; (b) the pressure amplitude of the pulse and the total energy fluence in the pulse, 
(pc)^^ dt. These last two quantities are plotted as functions of observation angle as defined in Figure ^. 



are presented in Figure ^. The multi-peak structure is 
clearly visible in the amplitude and fluence diagrams as 
function of the angle. While the pressure pulse shape, 
taken at the direction of maximum emission, has, also in 
this case, a simple bipolar shape, the maximum occurs 
at a large angle respect to the previous case, because of 
the longer shower profile. 

The signal from electron (anti)ncutrinos is the super- 
position of these two cases in a proportion correspond- 
ing to the way the energy of the primary is shared be- 
tween the hadronic and electromagnetic shower. It was 
found ||l9| that at £'^(p) — 10^° eV the hadronic compo- 
nent accounts for ^20% of the total energy. For an UHE 



primary, the axes of the hadronic and electromagnetic 
showers are practically parallel, the angle between them 
being of the order of 10~^ radians |34[ . 

We note here that the peak pressures predicted at 
1000 m from 10^° eV showers are well within the sensitiv- 
ity of good quality hydrophones. Indeed typical sensitivi- 
ties for the frequency band of interest are ~ 10"'^ Pa [ pO[ . 
However it is rather clear that the two factors limiting 
the power of this technique will be the ambient noise and 
the characteristic emission pattern described above. Such 
pattern substantially limits the solid angle accessible to 
each sensor and hence, together with the noise level, it 
will dictate the maximum tolerable spacing between de- 
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FIG. 5: The results of calculation of the acoustic signal from the electromagnetic part of a !/e(i7e)-induced shower at 
the distance of 1000 m from the shower axis. The total energy in the electromagnetic shower is 10'^° eV. (a) the pulse shape at 
the observation angle of 6' = 8°, where the amplitude is maximal; (b) the amplitude of the pulse and the total energy fluence 
in the acoustic pulse, (pc)^^ Both are plotted as functions of observation angle as defined in Figure ^. 
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III. THE AUTEC ARRAY AS A UHE 
NEUTRINO DETECTOR 

While the installation of hydrophones to complement 
underwater Cerenkov arrays has been discussed by sev- 
eral groups 1 45, it is interesting to consider whether 



a very large existing array could be used in parasitic 
mode for UHE neutrino detection. Although the scien- 
tific community has recently started discussing the use of 
long-range, early warnin g m ilitary arrays for a variety of 
oceanographic purposes , the low bandwidth of these 
systems, essentially designed to detect low-frequency ship 
noise that can propagate over very large distances, make 
them quite unsuitable for the type of signals discussed 
here. Another type of arrays exists, generally designed 
to track ships and weapons equipped with special high 
frequency "pingers" during limited-range naval exercises. 
Such arrays have a relatively high-density of hydrophones 
with typical bandwidths in the tens of kHz. The ar- 
ray considered here, at the Sites 3 and 4 of the Atlantic 
Undersea Test and Evaluation Center of the US Navy 
(AUTEC) j4|], covers an approximate area of 250 km^ 
(« 5 X 15 nautical miles^) with depths between 1400 and 
1600 m, as shown in Figure ^. 

AUTEC is located in the "Tongue of the Ocean", a 
50 X 200 km^ tract of deep sea, bounded to the west by 
Andros Island and to the south and east by large areas 
of very shallow banks in the Bahamas. This peculiar ge- 
ographical configuration, with shipping access only from 
the north through the narrow Providence Channel, pro- 
vides quiet conditions, because of the low boat traffic and 



sluggish currents. 

Individual hydrophones at Sites 3 and 4 are mounted 
at the ends of 4.5 m long booms extending above the 
ocean floor. A total of 52 sensors cover the two sites, 
arranged on a triangular lattice with 2.5 km sides. Ana- 
log signals from each hydrophone are preamplified and 
brought to shore where digitizers and processors are lo- 
cated. The frequency response of the hydrophone and 
analog chain is flat to within ±5 dB in the range from 
1 to 50 kHz, while the sampling rate of the digitizers is 
about 100 kHz. Accurate GPS time stamping is provided 
in the data stream. In its normal operation the system is 
capable to gather highly accurate 3-dimensional in- water 
tracking data. While a denser sensor spacing would of 
course be desirable, the shower-to-hydrophone distance 
of 1 km we use though this work is close to the worst 
case scenario of a neutrino interacting half-way between 
two sensors. 

In Figure ^ we show the random noise levels at AUTEC 
for different wind conditions, along with the cumulative 
probability for such conditions to occur on range. The 
approximate frequency spectrum of the expected neu- 
trino signals is also given for reference. The noise spectra 
in the Figure are due to the waves at the ocean surface 
induced by wind j4^, ^ ^ and, above >10 kHz, to ther- 
mal noise p9[ . In addition coherent noise from human 
activities and natural phenomena should be considered. 
A-priori we expect two types of artificial noise: the first, 
due to ship screws, is mainly confined to low frequency 
and hence easy to reject, while the second, due to the 
range "pingers" , is concentrated at a few well known fre- 
quencies and hence easy to filter-out. In general exer- 
cises are performed on range roughly 50% of the time, 
so that very little man-made noise is expected for a sub- 
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FIG. 6: Schematic view of the AUTEC facilities. The hy- 
drophone array discussed in the text is the larger quadrangu- 
lar region. 



stantial fraction of every day (typically outside of work- 
ing hours). More serious is probably the high frequency 
noise produced by marine mammals using sonar to local- 
ize their prey and snapping shrimps . The severity of 
these backgrounds depends on the season and can only 
be quantitatively understood by analyzing a substantial 
amount of data from the array. Here we limit ourselves to 
a detailed analysis of the random noise, formulating a fil- 
tering algorithm that can be implemented on the digital 
signal processors (DSP) that analyze on-line the time- 
series from each hydrophone. This system of data selec- 
tion can be used essentially as a trigger signal to log an 
interval of the data stream for the entire array (or maybe 
only some subset of sensors near to the one producing the 
trigger). Further data reduction, possibly involving the 
correlation of signals from different sensors, can then be 
done off-line. 

Scattering of sound off the ocean surface and bottom 
may alter somewhat the simple picture given above. For 
instance we expect that each hydrophone will generally 
record two pulses: the direct one and the one scattered 
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FIG. 7: Ambient noise spectra at AUTEC for different wind 
speeds (in knots, 1 knot = 0.5144 m/s) and the approximate 
spectrum of the expected neutrino signal. The absolute am- 
plitude of the signal is, of course, dependent upon the energy 
of the primary and the location of the impact relative to the 
hydrophone. The percent figure given next to each curve rep- 
resents the cumulative probability of finding such (or better) 
conditions at a given time. 



from the ocean floor immediately around the sensor. The 
time delay will be of the order ol h/c 3 ms, where 
h — 4.5 m is the height of the hydrophones above the 
bottom. This effect can be used to estimate the an- 
gle of incidence of the acoustic wave and hence of the 
shower axis. In addition scattering phenomena can make 
the event detectable by more than one hydrophone. The 
attenuation and re-emission patterns from scattering at 
the sea-bed and surface (where bubbles play an impor- 
tant role) can be numerically estimated ||5^, While 
large attenuations are to be expected so that it appears 
unlikely that these phenomena can be used to trigger an 
event that for geometrical reasons was not directly visi- 
ble, it is probable that scattered signals can be detected 
below threshold by sensors in the vicinity of the trigger- 
ing hydrophone. While off-line study of these correlations 
would be particularly important in confirming the event 
and in measuring position, orientation and energy of the 
shower, in the rest of this paper we concentrate on the 
triggering function that has to rely on single sensors. 



IV. THE SENSITIVITY OF THE DETECTOR 

In order to analyze the problem of signal detection 
in some detail we assume a discrete data sample set Xk 
(fc = 0, . . . , A/ — 1) with a sampling rate of fs = 100 kHz, 
as provided by the AUTEC digitizers. We consider the 
detection from a single site, as it is natural given the 
small probability that the radiation pattern discussed in 
the previous section would intercept more than one hy- 
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FIG. 8: The signal from 2 x 10^^ eV hadronic shower arti- 
ficiaUy embedded in the noise sample of length AI — 10000 
at position kg — 5000. The noise is simulated using the spec- 
trum for wind speed of 13 knots. The time interval of this 
sample is 0.1 s. 



drophone. Sub-threshold use of signals from neighboring 
hydrophones may be possible in off-line analysis but it 
is not discussed here. The data are the sum of the sta- 
tionary Gaussian noise of a given power spectrum and a 
signal of a given shape, but unknown amplitude, starting 
at position ka in the data set, as illustrated in the sim- 
ulated time series of length M = 10000 in Figure In 
this Figure a signal from a 2 x 10"'^^ eV hadronic shower at 
1 km with the optimal angular orientation to illuminate 
a hydrophone has been superimposed at a random time 
(chosen in this case to start at sample number ko = 5000) 
over the statistical noise spectrum relative to 13 knots 
wind from Figure ^ These conditions (or better) occur 
at the AUTEC site 60% of the time. The sampling fre- 
quency is high enough to provide negligible distortions. 

We extract the signal using the algorithm described in 
Appendix, which is based on the digital filter of transfer 
function Hi calculated on the basis of the signal shape 
and the Gaussian noise spectrum given in equation ( [A4| ). 
The variable Yfc resulting from the application of the filter 

k+N-l 

Yk = ^ HiXk+i 

l=k 

has a Gaussian distribution and is used to assert the pres- 
ence of the signal, using a threshold Yth that can be cho- 
sen for a certain detection efficiency and false-alarm rate 
at a given signal amplitude A {Yk > Ith indicating the 
presence of a signal). In practice the value of N can be 
made small enough for efficient calculation on the com- 
mercial DSP processors used at AUTEC. 

In Figure |^ we plot the variable Yfe for the time se- 
ries in Figure ^. We see that at the time sample 
k = ko = 5000 has a value substantially higher that 
elsewhere. Using a threshold Ith = Ykg, we find the 
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FIG. 9: The variable Y discussed in the text and computed 
in detail in Appendix (see equation A3) for the time series 
shown in Figure ^. A prominent peak at the correct location 
of fco — 5000 is clearly visible. 



probability of false alarm to be 5 x 10^^^ at each data 
point. On the other hand the naive technique of trying 
to find the signal using an amplitude threshold in time 
domain would give a false alarm probability of 0.14 at 
each sample point, with a threshold set at amplitude of 
the signal. The effectiveness of the method, also evident 
by simple inspection of Figures || and ^, is based on the 
difference between the signal and noise spectra in Fig- 
ure 1^. This technique is equivalent to the application of 



a matched filter M, 55 . A value of iV = 7 was used in 
the calculation, giving ~15 floating-point operations at 
each data point, which corresponds to a processor speed 
requirement of 1.5 Mflop/s. The duration and amplitude 
of the signal can be estimated by maximizing a likelihood 
function. 

We now use the probabilities of a signal miss and a false 
alarms found in the Appendix to quantitatively analyze 
the power of the array as a UHE neutrino detector. Here 
we stress the fact that our analysis does not take into 
account possible interference from coherent sources, as 
already mentioned above. 

We define the detector efficiency as the fraction of neu- 
trinos that are detected respect to those interacting in the 
volume of water £ = A/Aq. We can write the interaction 
rate Aq as 

Ao = F{E)NAPro'j{E)V^ 

where F{E) is the UHE neutrino flux at energy E, (j{E) 
is the cross-section of interaction with a nucleon, NaPw 
is the number of nucleons in unit volume of water, and 
Vw — ah is the volume of the detector (a and h being 
the area and average depth of the array). The detected 
rate is limited by the effective volume of the radiation 
pattern Ve in which the amplitude of the signal (relative 
to produced by a shower of energy 10^° eV at a distance 
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FIG. 10: False alarm rates A for different detection effi- 
ciencies £. The solid lines are for or (producing an 
electromagnetic and an hadronic showers) while the dashed 
lines are for other neutrino flavors (where only the narrower 
beam from an hadronic shower is detected). The energy in 
abscissa is relative to the neutrino and it is shared between 
hadronic and leptonic components as expected for the energy 
of 10^° eV (20% hadronic and 80% leptonic). 



of 1 km) is above the detection threshold A. We can 
write Ve as: 



Ve 
1 km^ 



27rA6' f Re 
3 

2nA9 



1 km 

E 



102" cV 



1 3 



(3) 



where A9 is the effective angular extent of the radiation 
lobe that can be found by integrating the curves in Fig- 
ures ||d and ||b, and Re is its radial extent. We obtain 
A9 ~ 0.5° {Ad ~ 1.5°) for showers induced by i^^, 9^, z/r, 
Vr (fe, Vf,). The difference comes about from the electron 
neutrinos having the electromagnetic part of the shower, 
while other flavors have only hadronic part. We note 
that the expression for Ve includes a signal amplitude 
inversely proportional to the distance i?o, as it should in 
our regime of radiation . 

The rate of detected neutrinos is then found by inte- 
grating over volume: 



A = NhF{E)Nap^<j{E) 



Ve,„ 



Pdetect{A,Yth) dVE 



(4) 

where Nh is the number of hydr oph ones and pdetect is the 
detection probability given in (A7), with A found from 
Ve using (^). The upper limit of the integration VE.max is 
determined from the condition that the acoustic radiation 
pattern is limited in size by the size of the detector, i.e. 



Re < -Rmax, where R„ 
V„ 



100 km. Finally, we get 



Pdetect{A,Yth) (We 



We can now use this expression to set the value of 1th 
for any given efhciency £ and then calculate the result- 
ing false alarm rate as A = NhfsPiedsc, where pfaise is 
determined from Yth using (A8). We plot the calcu- 



lated false alarm rates for different efficiencies in Fig- 
ure p^a. The solid lines correspond to or that pro- 
duce an electromagnetic and an hadronic showers, while 
the dashed lines are for other neutrino flavors that pro- 
duce a narrower sound lobe from the shorter hadronic 
shower. As expected electromagnetic showers can afford 
a higher threshold that results in a lower false alarm rate. 



V. CONCLUSIONS 

We have reviewed the possibility of searching for ultra- 
high-energy neutrinos in cosmic rays using the acoustic 
emission from the electromagnetic and hadronic show- 
ers produced by the neutrino interactions in sea water. 
We have analyzed the expected statistical noise and de- 
vised an optimal algorithm to filter it out of the data 
stream in real-time. In a simulation we have applied this 
technique to trigger the data-acquisition of a very large, 
high-frequency multi-hydrophone array of the US Navy 
and found that one could trigger on events at or above 
~ 10^° eV with tolerable false alarm rates. Our algo- 
rithms are based on signals from individual sensors, as 
appropriate for a trigger system, and are optimized to run 
on the array's digital signal processors. This would make 
a data taking campaign rather straightforward. Further 
off-line analysis would be needed to study the additional 
information that can be obtained by multi-hydrophone 
correlations. Although it is clear that an array with 
higher sensor density would allow lower energy thresh- 
olds and better redundancy, a test on the existent array 
analyzed here will provide information on coherent noise 
and give a definite assessment on the power of this tech- 
nique. 
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APPENDIX A: THE DETECTION ALGORITHM 

In order to analyze the problem of signal detection in 
some detail we assume a discretized data sample from 
a single hydrophone. Since the amplitude of the signal 
can vary with respect to the distance to the source, we 
only fix the shape of the signal, without specifying its 
amplitude in advance. Let us assume that the signal 
shape after discretization with sampling frequency /„, is 
represented as an array of N real numbers, Fn, where 
n = 0, . . . , iV — 1, so that the signal is AF^ where A is 
the amplitude. 

The discrete data sample Xk, fc = 0, . . . , Af — 1, which 
is discretized with sampling frequency /s, is the sum 
of the stationary Gaussian noise Wk of a given power 
spectrum and a signal embedded at a random position 

= . . . ,M ~ N with an unknown amplitude A: 



Xk = 



Wk 



k < ko OT k > ko + N 
ko < k < ko + N. 



We will search for the signal position fco and the am- 
plitude A which maximizes the conditional probability of 
a signal to be present within a given data set. This is a 
regression problem that can be solved by finding the max- 
imum in the appropriate likelihood function [p6[. Let us 
take a subsample of the original data Xk = Xko+k, where 
/c = 0, . . . , TV — 1, and analogously define a noise subsam- 
ple Wk = Wko+k- The likelihood of the signal presence at 
position ko is given by the Bayes formula 



Prob{F, Ae[A,A + dA]\X} 

_ Proh {X\F, A} PToh{F} f{A)dA 
^ PmhJX} 



(Al) 



The ratio to the probability of the absence of signal 
(which is found in a similar manner) is 

Prob {F, Ae[A,A + dA]\X} 
Probjno signaljX} 

Prob {X\F, A} Prob{F} fiA)dA 



Prob{X|no signal} Prob {no signal} 

We notice that some of these conditional probabilities are 
related to the distribution of the Gaussian noise: 



Proh{X\F,A} 
Prob{X|no signal} 



fw(X-AF)dX 
fw{X)dX 



where dX = Ilfc'^^fc- The probability distribution 
of the Gaussian noise, defined so that jw^^dW 
is the probability that the noise values Wk are in the in- 
terval dW = OfcLV ^W^fc' 

is given by 



jw{W) = Cexp ( - 2^ 



where K is the noise covariance matrix defined as Kjk = 
(VFjM^fc) that, for a stationary noise, has the property 
of being a function of the difference (fc — j) only. If this 
portion of the time series is long enough so that the noise 
values separated by > N (2 points are not correlated, 
we can assume circular stationarity Kjk — R{k-j) mod Nj 
where R is an array of length N. 

The calculations are easiest in a basis in which the 
covariance matrix is diagonal. For a circularly stationary 
signal, this is the basis of Fourier components 



E 

n=0 



Wk= ^ Wn exp 



2m, 
'1^ 



-kn 



0,...,7V-1 (A2) 



The distribution of the noise Fourier components is de- 
fined so that f-^{W)dW is the probability of W to be in 

the interval dW = dWodWM/2 nfi^^ d[ReWk]d[luiWk]. 
Hence the probability distribution is given by 



/,;^(W) = C'exp(^-i^|V7fc|V5fej 

where Sk, k = . . . ,N — \ is the noise spectrum defined 
as {W^Wi) — SkiSk- This is related to the covariance 
matrix as 5 = NR, with R defined above. 

Substituting f^^{W) we get the ratio of the likelihoods 



Prob {F, Ae[A,A + dA] \X} 

Probjno signaljX} Prob {no signal} 

where the argument of the exponential function is 



\Xk-AFk\^-\Xk\^ 



In /(A) 



Instead of maximizing the original probability, we maxi- 
mize the log-likelihood L' . 

Since A is always positive and \nf{A) w const for A > 
0, wc choose the optimal amplitude estimate A which 
maximizes L': 



A = max < 0, 



EkMXkF^)/Sk 
EklhWSk 



Substituting this value into L' , we form a random vari- 
able 

L = -IJ2^ (l^fc - - \Xk\^) - L'{A) + const. 

k 

Note that the maximization with respect to A can be 
generalized to other parameters of the signal. Let F*^™^ 
be different signal shapes from a set given by an index m. 
Then to get the maximum likelihood of the signal shape 
we choose the value of mg such that = max{L('")} 

(note that const = \nf{A) is the same for all m). 
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It turns out that the equation for L can be simplified 
in terms of the number of calculations required. Substi- 
tuting into it the expression for A, we derive 



L 



where Y = Efc Re(Xfcf^*)/5fe and al = Y^kW/Sk 
can be shown to be the vari ance of Y . Substituting the 
Fourier transform definition (A2) for X, wc finally obtain 




JV-l 



1=0 



where 



(A3) 



(A4) 



The variable Y is thus the result of applying a digi- 
tal filter with response function H or transfer function 
Hk = NFk/Sk- The amount of calculation can be re- 
duced even further, due to a fact that the quantity of 
Hk that are significantly different from zero is N* < N. 
The number TV* is estimated to be ~ Tgig/s because Hk 
has the time duration of the same order as the duration 
of the signal, Tsig. This last fact is confirmed by a di- 
rect calculation. Thus, for typical values of fs = 10^ Hz 
and Tsig — 10""* s, we have N* ~ 10. The number of 
floating-point operations (flops) required is then ~ 10^ 
per second, which is feasible on modern DSP processors 
capable of more that 50 Mflops/s. 

Signal detection can then be implemented at every hy- 
drophone utilizing the fact that Y has different distribu- 
tions fY,A{Y) in the cases presence {A ^ 0) and absence 
{A = 0) of a signal. We chose a threshold value Yth and 
decide that the signal is present (absent) when Y > Fth 
(Y < Yth)- Since the distributions fY,A{Y) for the cases 
of signal presence and absence in general overlap, there 
will be the possibility of false alarms and signal misses. 



The threshold value has then to be chosen in an appro- 
priate manner. 

We now estimate the probabilities of false alarms and 
misses for a certain threshold Yth by estimating the prob- 
ability distribution fY,A{Y) of the variable Y in the pres- 
ence of a signal of amplitude A (including the no-signal 
situation as a particular case ^ = 0). Since y is a lin- 
ear combination of Gaussian variables, the distribution 
sought is also Gaussian 



1 g-(y-Aa^)V(2<T^)^ 



The probability of missing a signal is then 



Pn 



fYAY)dY 

-oo 

- 1 - ^erfc r^^-^^- 
while the probability of the signal detection is 



(A6) 



Pdotoct 1 - Pn 



fYAy)dY 



-erfc 



1 „ fY,^-Aal 



2""" \ ^CTY 

Finally, the false alarm probability is 
Pfaisc = / fY,A=o{Y) dL = -eric 



V\/2, 



cry 



(A7) 



(A8) 



In the expressions above erfc (x) = 1 — erf (x) is the com- 
plementary error function. The above expressions can be 
inverted to derive Fth from pdctcct and pfaisc: 



_ / Aalr + V2cryerf"^ (1 - 2p 

detect ) 

" 1 V2ai'erf-i (1 - 2pfaisc) 
where erf~^ (x) is the inverse error function. 



(A9) 
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